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Abstract 

We perform a dispersive theoretical study of the reaction 77 — > 7r°7r° emphasizing the low energy 
region. The large source of theoretical uncertainty to calculate the 77 — > 7r°7r° total cross section for 
\fs > 0.5 GeV within the dispersive approach is removed. This is accomplished by taking one more 
subtraction in the dispersion relations, where the extra subtraction constant is fixed by considering 
new low energy constraints, one of them further refined by taking into consideration the /o(980) region. 
This allows us to make sharper predictions for the cross section for y/s < 0.8 GeV, below the onset 
of D-wave contributions. In this way, were new more precise data on 77 — > 7r°7r° available one might 
then distinguish between different parameterizations of the tttt isoscalar S-wave. We also elaborate on 
' the width of the a resonance to 77 and provide new values. 
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1 Introduction 



The reaction 77 — > 7r°7r° measured in ref.[lj offers the interesting prospects of having a two-body hadronic 
final state and the important role of final state interactions in S-wave enhanced due to the null charge 
of the 7T°. These two facts make this process very suited for learning about the non-trivial 7T7t isospin 
(I) S-wave. In addition, it was taken as an especially appropriate ground test for Chiral Perturbation 
Theory (xPT) [21 [3], since at lowest order this process is zero and at next-to-leading order (one loop) is 
a prediction free of any counterterm [H[5]. However, the one loop xPT prediction departs very rapidly 
from data just above the threshold and only the order of magnitude was rightly foreseen. A two loop 
calculation in ref . [U [7] was then undertaken with better agreement with data pQ . The three counterterms 
that appear at 0(p 6 ) are fixed by the resonance saturation hypothesis. Other approaches supplying higher 
orders to one loop xPT by taking into account unitarity and analyticity followed [El EJ [TO] . Ref. [TO] is a 
Unitary xPT calculation in production processes [TOl [HI [121 [131 [HI [15] and was able to provide a good 
simultaneous description of 77 — ► 7r°7r°, 7r + 7r~, mr°, K + K~ and K°K° from threshold up to rather high 
energies, s 1 / 2 < 1.5 GeV. This approach was also used in ref. [16] to study the 77 — > 7r°77 decay. 

We concentrate here on the dispersive method of refs. |17|. \18\ 119] , We critically review and extend 
it, so as we are able to drastically reduce the uncertainty due to the not-fixed phases above the KK 
threshold of the 1 = S-wave 77 — ► tttt amplitude. This is accomplished by using an I = S-wave 
Omnes function that is continuous under changes in the phase function employed for its evaluation 
above the KK threshold. Equivalently, one can introduce an additional subtraction in the dispersion 
relation to evaluate the 1 = S-wave 77 — ► tttt amplitude to those considered in ref. [El HHJ HHJ [8]. 
This new subtraction constant is fixed by considering simultaneously three constraints instead of the two 
employed in the previous references. As a result of this much reduced uncertainty, the total cross section 
<t(77 — ► 7r°7r°) might be used to distinguish between different S-wave 1 = phase shift parameterizations 
once new more precise experimental data become available. We also perform calculations of the width 
r(ci — > 77), taking from the literature different a resonance parameters, and compare with the value of 
ref.[T7j. Other papers dedicated to calculate the two photon decay widths of hadronic resonances are 

The content of the paper is as follows. In section [2] we discuss the dispersive method of ref.[18j and 
extend it to calculate with higher accuracy the cross section 77 — ► 7r°7r°. The resulting 17(77 — * 7r°7r°) 
and r(cr — * 77) are given in section [3l We elaborate our conclusions in section [H 

2 Dispersive approach to 77 — > 7r°7r° 

In refs.|17 | 118 1 H9j an interesting approach was established to calculate in terms of a dispersion relation 
the 77 — * (tttt)j S-wave amplitudes, Fj(s), where the two pions have definite isospin I. Notice that for 
77 —* 7r°7r°, due to the null charge of 7r°, there is no Born term, fig[TJ One then expects, as remarked 
in ref.[l8j, that only the S-wave would be the important partial wave at low energies, < 0.7 GeV. 
For 77 — ► 7r + 7r~, where there is a Born term due to the exchange of charged pions, the D-waves have 
a relevant contribution already at rather low energies due to the smallness of the pion mass. In the 
following, we shall restrict ourselves to the S-wave contribution to 77 — > 7r°7r°. The explicit calculation of 
ref. [10] indicates that the D-wave contribution at y/s ~ 0.65 GeV is smaller than a 10% in (1(77 — ► 7r 7r ), 
and it rapidly decreases for lower energies. 

The function Fj(s) is an analytic function on the complex s— plane except for two cuts along the real 
axis. The right hand cut happens for s > 4m 2 , with m n the pion mass, and is due to unitarity. The left 
hand cut, in turn, runs for s < and is due to unitarity in crossed channels. Let us denote by Lj(s) the 



2 





Figure 1: Born term contribution to 7(^:1)7(^2) — ► P + (pi)P~(p2)- 

complete left hand cut contribution. Then the function Fj(s) — Lj(s), by definition, only has right hand 
cut. Next, refs.fTT l \18 \ [T9] consider the Omnes function cjt(s), 



o;/(sJ = exp 



(2.1) 



with 4>i(s) the phase of Fj(s) modulo ir, chosen in such a way that (f>i(s) is continuous and 0/(4m^.) = 0. 
Because of the choice of the phase function ^j(s) in eq. (|2.1j) . the function Fi(s)/ujj(s) has no right hand 
cut. Then refs.|17 1 1181 119j perform a twice subtracted dispersion relation for (Fi(s) — Li(s))/u>i(s), 

)=L/(s +«/^/(« ) + cjsw 2 (*) + —"/ (*) / ,2! , m / ,M ds • 2 - 2 

7T 74 m 2 S M (s'- S)\U)1[S') I 

On the other hand, Low's theorem [21] requires Fj — > Bj(s) for s — ► 0, with 5/ the Born term contri- 
bution, shown in figHJ If we write Lj = Bj + it!/, with 72/ — > for s — * 0, as can always be done, then 
Low's theorem implies also that Fj — Lj — > for s — ► and hence ai = 0. 

For the exotic 1 = 2 S-wave one can invoke Watson's final state theorenj^l so that 4>2{s) = S n (s)2, the 
1 = 2 S-wave tttt phase shifts. For 1 = the same theorem guarantees that 4>o(s) = 6 w (s)o for s < Am K , 
with 5- k {s)q the S-wave 1 = irir phase shifts and itik the kaon mass. Here one neglects the inelasticity 
due to the 4-7T and 6ir states below the two kaon threshold, an accurate assumption as indicated by 
experiment |22[ 123] . Above the two kaon threshold, sk = 4m 2 ^, the phase function 4>o(s) cannot be fixed 
a priori due to the onset of inelasticity. This is why ref.[l7] took for s > sk that either (j>o(s) — S n (s)o 
or (f>o(s) ^ 5 n (s)o — 7T, in order to study the size of the uncertainty induced for low energies. It results, 
however, that this uncertainty increases dramatically with energy such that for = 0.5, 0.55, 0.6 and 
0.65 GeV it is 20, 45, 92 and 200%, see fig.3 of ref.[H]. 

The reason for this behaviour is the use of the function ujo(s) in eq. (|2.2p . The 1 = S-wave phase 
shift <5ti-(s)o has a rapid increase by n around sk = 4771^, due to the narrowness of the /o(980) resonance 
on top of the KK threshold. Let us denote by <p(s) the phase of the irir — * irir 1 = S-wave strong 
amplitude, modulo ir, such that it is continuous and tp(4m 2 ) = 0. This phase is shown in figE] together 
with 5jr(s)o and ^(5)2- Now, if one uses (p(s) instead of (f>o(s) in eq. (|2.ip for illustration, the function 
u>o(s) is discontinuous in the transition from 5k(sk)o — > vr — e to 5 n (sK)o — ¥ ^r + e, with e — ► + . In the first 
case |wo(s)| has a zero at sk, while in the latter it becomes +00. This discontinuity is illustrated in fig[3] 
by considering the difference between the dot-dashed and dashed lines. This discontinuous behaviour of 
loo(s) under small (even tiny) changes of 5 n (s)o around sk, was the reason for the controversy regarding 
the value of the pion scalar radius (r 2 )^ between [24 1 125 1 126] and [27]. This controversy was finally solved 



#1 This theorem implies that the phase of Fi(s) where there is no inelasticity is the same, modulo ir, as the one of the 
isospin I S-wave tttt elastic strong amplitude. 
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in ref.[28j where it is shown that Yndurain's method is compatible with the solutions obtained by solving 
the Muskhelishvili-Omnes equations for the scalar form factor |29 |, 1301 f3Tj . The problem arose because 
refs.[2U ES] overlooked the proper solution and stuck to an unstable one. 




Figure 2: The phase shifts 6V(s)o and 6V(s)2 and the phase ip(s). Experimental data are from refs. [3"2l [251 [551 IM] 
for 1 = and refs435, 36J for 1 = 2. The insert is the comparison of CGL [37] and PY [38] with the accurate data 
from K e 4 [331 OS]. 
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Figure 3: |o>o(s)|, eq . (|2. ljl . with 6V(s_k-)o < tt, dashed-line, and iv(sif)o > t, dot-dashed line. The solid line is 
|fio(s)|, eq. (j2.3| . for the latter case. Here ip(s) is used as ^o(s) in eqs. (j2.ip and (j2.3p for illustrative purposes. 

Inelasticity is again small for 1.1 < ^/s < 1.5 GeV being compatible with zero experimentally [22, 23J. 
As remarked in refs. [241 128 j . one can then apply approximately Watson's final state theorem and for 
Fq(s) this implies that 4>o(s) ~ <5W(s) modulo tt. Here d^ + \s) is the eigenphase of the irn, KK 1 = 
S-wave S-matrix such that it is continuous and S^ + \sk) = ^^(s^o- I n refs. [25] [28] it is shown that 

o - 7r, depending on whether (5 vr (s_ft:)o > vr or < n, respectively. In order to 
fix the integer factor in front of tt in (f>o(s) ~ <5W(s) modulo n, one needs to devise an argument to 
follow the possible trajectories of </>o(s) in the narrow region 1 < ^/s < 1.1 GeV, where inelasticity is not 
negligible. The remarkable physical effects happening there are the appearance of the /o(980) resonance 
on top of the KK threshold and the cusp effect of the latter that induces a discontinuity at sk in the 
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derivative of observables, this is clearly visible in figi2J Between 1.05 to 1.1 GeV there are no further 
narrow structures and observables evolve smoothly. Approximately half of the region between 0.95 and 
1.05 GeV is elastic and tfio(s) = 5 n (s)o (Watson's theorem), so that it raises rapidly. Above 2rriK — 1 GeV 
up to 1.05 GeV the function 4>q(s) can keep increasing with energy, as 5 n (s)o or tp(s) for S^(sk)o > ^, 
and this is also always the case for the corresponding phase function of the strange scalar form factor 
of the pion [28]. It is also the behaviour for 4>o(s) corresponding to the explicit calculation of ref.|10j. 
The other possibility is a change of sign in the slope at sk due to the KK cusp effect such that 4>o(s) 
starts a rapid decrease in energy, like ip(s) for 5 w (sk)o < figGJ Above yfs = 1.05 GeV, <f>o(s) matches 
smoothly with the behaviour for yfs > 1.1 GeV where it is constraint by Watson's final state theorem. 
As a result of this matching, for > 1 GeV either 4>q(s) ~ 5 n (s)o or c/>q(s) ~ 5 n (s)o — ir, corresponding 
to an increasing or decreasing <fto(s) above sk, respectively. There is then left an ambiguity of ir in <f)o(s) 
for 1.5 GeV > y/i > ^/s~k- Our argument also justifies the similar choice of phases in ref.[T7] above sk 
to estimate uncertainties. Let us define the switch z to characterize the behaviour of 4>o(s) for s > sk 
such that z = +1 if 4>o(s) rises with energy and z = —1 if it decreases. Above 1.5 GeV the phase function 
employed has little effect in our energy region and we use the same asymptotic phase function as in 
ref.[28], tending either to 2tt (z = +1) or tt (z = —1) for s — > +oo. It allows a large uncertainty of ~ 2tt 
at yfs = 1.5 GeV, that only shrinks logarithmically for higher energies. This uncertainty is included in 
our error analysis. Further details are given in ref. |28j . 

Next, we define, the function 0,q(s), similarly as done in ref.[28j, 

n («) = U - d ( z )j-) ex P 

where 9{z) = +1 for z = +1 and for z = —1 and s\ is the point at which </>o(si) = tt. The latter is the 
only point where the imaginary part of ^o( s ) vanishes around sk and this fixes the position of the zero. 
Now, we perform the same twice subtracted dispersion relation as in eq. (|2.2f) but for {Fq{s)—Lq{s))/Q.q{s) , 

F (s) = L (s) + c sn (s) + S -n (s) f° f° ( f } ^/fl ds> + 9(z)^f\^(F ( Sl ) - L ( Sl )) . (2.4) 

J Aral S ( S ~ S ) I °0 («') | U>0{S1)S( 

In the previous equation we introduce <Pq(s) that is defined as the phase of fio(s)- Let us note that in 
the case z = +1 the phase of f2o(s) f° r s > si is not 4>o(s) but 4>o(s) — tt, due to the factor 1 — (s + ie)/s\ 
in f2o(s), eq. ()2.3p . Since 4>2(s), because of Watson's final state theorem, is given by S n (s)2, which is small 
and smooth [33 ES]) figEl the issue of the discontinuity in u>2(s) under changes in parameterizations 
does not rise and we use the dispersion relation in eq, (|2.2p . It is worth mentioning that our eq. (|2,4p for 
z = +1 is equivalent to take a three times subtracted dispersion relation for (Fq(s) — Lq(s))/loo(s), two 
subtractions are taken at s = and another one at s\. We could have taken the three subtractions at 
s = 0, although we find more convenient eq. (|2.4H which is physically motivated by the use of the Omnes 
function eq. (|2,3p that is continuous under changes in the parameterization of the 1 = S-wave S-matrix. 
When eq. (|2.3p is used with (p(s) instead of (fro(s) for 5 n (sK)o > 7r the solid curve in figj3] is obtained, 
which is again close to the dashed line for 5 n (sK)o < v. 

We denote by Fj^(s) the S-wave 77 — ► 7r 7r° amplitude and by Fc(s) the 77 — > tt + it~ one. The 
relations between Fq, F2 and Fn(s), Fc(s) in our isospin convention are 

F N (s) = -^=F + y|F 2 , F c (s) = -^=F - ^F 2 . (2.5) 

We have the unknown constants Co, C2 and Fo(si) — Lq(s\), the latter for z = +1. To determine them 
we impose: 



i r Ms') ds , 

7T 7 4m 2 S'(s' - S) 
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1. Fc(s) — Bc(s) vanishes linearly in s for s — * and we match the coefficient to the one loop %PT 
result [HE]. 

2. Fjsf(s) vanishes linearly for s — > as well and the coefficient can be obtained again by matching with 
one loop xPT 01(5]. 

3. For 1 = and z = +1 one has still Fq(si) — Lq(s\). The value of this constant can be restricted taking 
into account that -F/v(s) has an Adler zero, due to chiral symmetry. This zero is located at sa = hi 
one loop %PT and moves to sa = 1.175 in two loop %PT [6]. This implies about a 20% correction, 
that prevents us from taking a definite value for sa- In turn, we obtain that the value of the resulting 
cross section 17(77 ~~ > 

) around the /o(980) resonance is quite sensitive to the position of the Adler 
zero, because it controls the size of Fq(si) — Lq(si). The latter appears in the last term in eq. (|2.4p . the 
one that dominates Fq(s) around the /o(980) position since fio(si) = 0. Though the dispersive method 
is devised at its best for lower energies, it is also clear that it should give at least the proper order of 



magnitude for (1(77 —* vr vr ) in the /o(980) region !^ I Being conservative, we shall restrict the values 



of Fq(si) — Lq(si) so as the cross section at s\ is less than 10 times the experimental value around the 
/ (980) region, £7(77 -» 7r°7r°) < 400 nb. 

Regarding Lj(s), it is expected to be dominated at low energies by the Born term in isospin / because 
of the smallness of the pion mass. The Born term originates by the exchanges in the t and u channels 
(77r — * 77r) of charged pions, figHJ Other crossed exchanges of vector and axial vector resonances are 
relatively suppressed for Fj(s) due to the larger masses of the members of the J p = 1 , 1 ++ and l" 1 
multiplets so that their associated left hand cut singularities are further away. Among them, the 1 ++ 
axial vector exchange contributions are the dominant ones in 771^ — ► 'yir^ and already appear at the one 
loop level in xPT- The 1 and l" 1 exchanges start one order higher. As already remarked in ref.[8], 
the authors of refs. |18|, I19j. and later on also ref.[T7], overlooked the axial vector exchange contributions 
altogether and hence they are missing an essential part in the study of the low energy 77 — ► 7r°7r°. Indeed, 
once the 1 ++ axial vector exchange contributions are considered the cross section in the latter references 
would increase significantly at low energies. For instance, at y/s = 0.5 GeV one has more than a 20% 
increase compared to the case in which only the Born term and the 1 vector resonances exchanges 
are considered. At low energies the influence of the l" 1 axial vector nonet is much smaller than that 
of the 1 ++ and 1 multiplets. In the following, L[(s) is modelled by the Born terms and the crossed 
exchanges of the 1 ++ , 1 , and l" 1 resonance multiplets, evaluated from chiral Lagrangians. Explicit 
expressions for the distinct contributions to Lj(s) will be given elsewhere j 1 T| . 



3 Results 

In this section we show the results that follow by the use of eq. ()2.4p for 1 = and eq. (12.2p for 1 = 2. 
Since the main contribution in the low energy region to f2o( s ) an d ui2(s) comes from the low energy 
7T7T phase shifts, one needs to be as precise as possible for low energy irir scattering data. The small 
1 = 2 S-wave irir phase shifts, which induce small final state interaction corrections anyhow, can be 
parameterized in simple terms and our fit compared to data can be seen in figEJ For the 1 = S-wave 
7T7T, we take the parameterizations of ref.[37] (CGL) and ref.[38] (PY). Both agree with data from K e 4 
decays [331 El] and span to a large extend the band of theoretical uncertainties in the 1 = S-wave -kit 
phase shifts [ZQ [23102]. PY, similarly to refs. [4*2l |4"3"1 14"4"]. runs through the higher values of 5^(s)o, while 
CGL does through lower values, see figlH We shall use CGL up to 0.8 GeV, since this is the upper limit 
of its analysis, and the K-matrix of Hyams et al. [22\ above that energy. The latter corresponds to the 
energy dependent analysis of the experimental data of the same reference. On the other hand, PY is 

#2 E.g., other models |10l 1391 HO] having similar physical mechanisms describe that energy region very well indeed. 
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used up to 0.9 GeV, since at that energy this parameterization agrees well inside errors with [22], and 
above 0.9 GeV the K-matrix of ref.[22j is taken. Given the input functions (j>i(s) and Lj(s), the constants 
C2, Co and Fq(si) — Lq(si) can be fixed by the three conditions explained at the end of section [2l The 
77 — > 7T°7r S-wave amplitude Fn(s), eq. H2.5j) . can then be calculated and the total cross section is given 
7r°7r°) 



by cr(77 



with P(s) = \/l — 4m* /s. We show in figJUthe drastic reduction in 



the uncertainty of the cross section due to the variation of (j>o(s) above sk as commented in the previous 
section. For z = +1 one has the solid line while for z = —1 the dashed line results, both are very close. 
This should be compared with the dot-dashed line that is obtained from the approach of refs.|17 j [T8 | ll9|. 
The uncertainty now, by employing eq. (|2.4|) instead of eq. (|2.2p for I = 0, is drastically reduced. This 
improvement also implies that our results can be compared with data for y/s > 0.5 GeV. We also show 
by the gray band around the solid line the mild influence in our calculations of the uncertainty in the 
location of the Adler zero, restricted so that 17(77 — * Tr°n°) < 400 nb at the /o(980) region (experiment 
is ~ 40 nb). 
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Figure 4: The solid line corresponds to z = +1 and the error band is the uncertainty by requiring that (7(77 — * 
7r°7r°) < 400 nb at s±. This line should be compared with the dot-dashed one that would result from the formalism 
of ref.|17j. including axial vector exchanges. Finally, the dashed line corresponds to z = — 1. 

Our final (7(77 — ► vr 7r ) is shown in figJSJ We give the corresponding results for CGL (solid) and 
PY(dashed), where the band around every line stems from the uncertainties in our approach, which 
comprise: the errors in the Hyams et al. [22], CGL and PY parameterizations (those of the last two 
indeed dominate the width of the bands), to use either 0q(s) ^ <5,i-(s)o or 5 n (s)o — n for s > sk, the 
uncertainty in the asymptotic phase and to restrict (1(77 — ► 7r°7r°) < 400 nb in the /o(980) region for 
z = +1. On top of that, we evaluate the conditions 1 and 2 above from the expressions given by one loop 
xPT either by employing f n = 92.4 MeV or / ~ 0.94/„-, where the former is the pion decay constant and 
the latter is the same but in the SU(2) chiral limit [3]. This amounts to around a 12% of uncertainty in 
the evaluation of Co and C2, due to the square dependence on /„.. Note that both choices, f n or /, are 
consistent with the precision of the one loop calculation and the variation in the results is an estimate for 
higher order corrections. However, the error induced in (7(77 — ► 7r°7r°) is much smaller than that from 
the other sources of uncertainty and can be neglected when added in quadrature. 
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In the same figure the dotted line corresponds to one loop CHPT [H |5] and the dot-dashed one to the 
two loop result [6j [T] . The latter is closer to our results but still one observes that the 0(p 8 ) corrections 
would be sizable. It is worth stressing that if the axial vector exchanges were removed, as in refs. |17|, 119] . 
then our curves would be smaller. This corresponds to the dot-dot-dashed line in figE] which is very close 
to that of ref.[T7] when employing c/>q(s) ~ S w (s)o — ir for s > sk- This curve is evaluated making use of 
CGL and ref.[22j. The three experimental points [T] in the region 0.45 — 0.6 GeV agree well with this 
curve. However, once the axial vector are included the curve rises. These three points lie around 1.5 
sigmas below the CGL result band, and by more than two sigmas below the PY one. This clearly shows 
that more precise experimental data on 77 — ► 7r°7r could be used to distinguish between different S- 
wave parameterizations. In turn, the next three experimental 17(77 ~~ > ^ ^°) points, those lying between 
0.6 — 0.75 GeV, agree better when the axial vector resonance contributions are taken into account, as one 
should do. As a result of this discussion, more precise experimental data for 77 — ► 7r°7r° are called for. 








Figure 5: Final results for the 77 — ► 7r°7r° cross section. Experimental data are from the Crystal Ball Coll. pQ, 
scaled by 1/0.8, as | cos#| < 0.8 is measured and S-wave dominates. The solid line corresponds to CGL and the 
dashed one to PY. The dot-dot-dashed line results after removing the axial vector exchange contributions. The 
band along each line represents the theoretical uncertainty. The dotted line is the one loop %PT result [?J [5] and 
the dot-dashed one the two loop calculation [BJ. 

In terms of the calculated Fn(s) one can evaluate the a coupling to 77, called g CT77 . The dispersion 
relation to calculate Fq(s) is only valid on the first Riemann sheet. If evaluated on the second Riemann 
sheet there would be an extra term due to the a pole. However, the relation between Fq(s) and Fq(s), 
the latter on the second sheet, can be easily established by using unitarity above the tttt threshold, 

F (s + ie) -F (s- ie) = -2iF (s + ie)p(s + ie)T?j(s - ie) , (3.6) 

with 4to^ < s < 4m|-, p(s) = (3(s)/1Qtt and e — > + . In the equation above Tjj(s) is the 1 = S-wave tttt 
elastic amplitude on the second Riemann sheet and Tj(s) is the one on the physical Riemann sheet. Due to 
continuity when changing from one sheet to the other, Fq(s— ie) = Fo(s+ie) , Tj =0 (s—ie) = Tjf Q {s+ie) . 
Then, eq. (|3.6p can be rewritten as 

F (s) = F (s)(l + 2ip(s)Tjf (s)) . (3.7) 
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Around the a pole, s a , 



T i= = ; Fo(s) = ^ ^7fa ^ (3 8) 



with QfjjfjY the <T coupling to two pions such that T = \g a7T7V \ 2 f3 /16irM , for a narrow enough scalar resonance 
of mass M. Notice as well the ^/2 factor in Fq(s) to match with the g^ir normalization used (the so 
called unitary normalization \A2 \ Wi \ H4]). Then from eqs. (|3.7|) and (|3.8|) it follows that 

2*77 _ 1^( Sct )V Fo(S(t)2) (3 _ 9) 



9l™ 2 V Stt 

Let us stress that this equation gives the ratio between the residua of the S-wave / = 77 — > tttt and 
7T7T — > tttt amplitudes at the a pole position. 

In order to derive specific numbers for the previous ratio in terms of our dispersive approach one needs 
to introduce s a . We take two different values for s a = (M a — iT a /2) 2 . From the studies of Unitary xPT 
[42| l43l IIIl 145] one has M a and T a around the interval 425-440 MeV. The other values that we will use 
are from ref. |46j . M£ cl = 441^g 6 MeV and T^ 1 = 544^25 MeV, where the superscript ccl indicates, in 
the following, values that employ the a pole position of ref.[l6]. The corresponding ratios of the residua 
given in eq. (!3.9p are: 



9a 



77 



$0-77 



(2.10 ± 0.25) x 10" 3 , S/j from ref.[ 

(2.06 ± 0.14) x 10 -3 , s a from ref.06] . (3.10) 

1 /2 

Both numbers are very similar despite that the imaginary parts of the two sj differ by ~ 20%. The 
result of |17j . with which we shall compare our results later, corresponds to the ratio in eq. (l3.10p being 
20% bigger at (2.53 ± 0.09) x lO" 3 with s a of ref. 06]. 

These ratios of residua at the a pole position are the well defined predictions that follow from our 
improved dispersive treatment of 77 — ► (tttt)i. However, the radiative width to 77 for a wide resonance 
like the a, though more intuitive, has experimental determinations that are parameterization dependent. 
This is due to the non-trivial interplay between background and the broad resonant signal. An unam- 
biguous definition is then required |17} 119] . We employ, as in ref. [17], the standard narrow resonance 
width formula in terms of g a ~ n determined from eq. (|3.9p by calculating the residue at s a , 

Nevertheless, the determinations of the radiative widths from this expression and those from common 
experimental analyses can differ substantially. The following example makes this point clear. 

From ref.[l5] one obtains \g a7r7T \ = 2.97 — 3.01 GeV, corresponding to the square root of the residua of 
the 1 = S-wave tttt amplitude, as in eq. (|3.8p . If similarly to eq. (|3.1ip . one uses the formula, 



the resulting width lies in the range 309 — 319 MeV, that is around a 30% smaller than T a ~ 430 MeV 
from the pole position of ref.[l5]. This is due to the large width of the a meson which makes the I^o-ttttI 
extracted from the residue of Tj^ , eq. (|3.8p . be smaller by around a 15% than the value needed in 
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eq. (|3.12p to obtain T CT ~ 430 MeV. Similar effects are then also expected in order to extract T(a — ► 77) 
from the eq. (|3.1ip . Equations similar to this are usually employed in phenomenological fits to data, e.g. 
see ref . |47j . but with |<7a- 77 | determined along the real axis. As a result of this discussion, one should allow 
a (20 — 30)% variation between the results obtained from eq, (|3.1ip and those from standard experimental 
analyses that still could deliver a 77 —* tttt amplitude in agreement with our more theoretical treatment 
for physical values of s. 

We shall employ the following values for \g a7rn \. First we take |# CT7r7r | = 2.97-3.01 GeV [4"2l I331 14"4" 1 14"5] . 
With this value the resulting two photon width from eqs. (|3.10p and (|3.1ip is 

r(<7 ->• 77) = (1.8 ± 0.4) KeV . (3.13) 

We also consider a larger value for |<7 CT7r7r | since T^ 1 [46] is larger by a factor ~ 1.3 than T a from ref. [45]. 
One value is 



Smrm - |fl™| "FT r = ( L127 ± °- 022 ) 1^1 = ( 3 - 35 ± °- 08 ) GeV ■ ( 3 - 14 ) 



T ccl (a -> tttt) 
(i) ~ \9°™\ \ T(a^7nr) 

This corresponds to the scenario discussed previously in eq. (|3.12| ) with a value 15% lower than 

1 fa,i d (2 ) = ^l^fj 1 ) ' = (3 - 93 ± °- 08) GeV ' (3 - 15) 

obtained by reproducing V£ l from the pole position using eq. (|3.12j) . If we evaluate with these couplings 
the <j — ► 77 width one obtains from eqs. (|3.10p and (|3.1ip . respectively, 

rg(«r->77) = (2.1 ± 0.3) KeV , 

rg(<r->77) = (3.0 ± 0.3) KeV . (3.16) 

Recently, ref.[T7] calculated a value T(a — > 77) = (4.09 ± 0.29) KeV also employing s CT from ref.[46j. 
This value is larger than T%Ua — ► 77) in the previous equation, despite that \g a7T7r \ there used is 3.86 GeV, 
very close to |<7<77r7r|(2y ^ is worth stressing that both our eq. (|3.1ip and eq.(7) of ref.[l7] are equivalent for 
calculating T(a — > 77). except for an extra factor |/3(so-)| ~ 0.95 in ref.[l7]. Of course, they are written 



in a different notation]^ 3 ! The reason for this remaining difference is two fold. As already mentioned 
above, ref. [17] does not include axial vector exchanges in evaluating 77 — * (tttt)i. It is this omission that 
accounts for half of the 20% difference in the ratio of residua, eq. (j3.10| ). mentioned above. The other 10% 
comes from improvements delivered by our extra subtraction and our slightly different inputs. Using the 
same value for \g^ n \ as in [17], our resulting value for T(a — ► 77) would be around a 40% smaller (as in 
eq. ipnfi]) ) than that in [17] . 

As a summary of the a — > 77 considerations, from our dispersive approach and s a of refs. |45l [35] 
we obtain a value for the ratio of the residua \g all /g anw \ ~ (2.1 ± 0.25) x 10 -3 . This number follows 
unambiguously from our study. Other more intuitive, but convention dependent quantities, like the 
<j — ► 77 width calculated from eq. (|3.1ip . are less well determined. These depend critically on the input 
value for |<7 CT7r7r | 2 and s CT , though they are not required in our dispersive study of 77 — ► 7r°7r°. We then 
determine the values: i) T(a — ► 77) = (1.8 ± 0.4) KeV with s a and Igawnl ~ 3 GeV from ref.[35]; ii) 
rg(cr -> 77) = (2.1 ± 0.3) KeV and Tp,(a -» 77) = (3.0 ± 0.3) KeV which come by considering s a 
of ref.[46j with an estimated \g a Tm\ — 3.4 and 3.9 GeV, respectively. Other values could be obtained 
from eq. ([3.11|) by plugging different s a and \g a nn\ in eq. (|3.9p in order to estimate |po- 77 |. One should 
require that these values are provided from a tttt S-wave 1 = strong amplitude in agreement with the 
experimental phase shifts, see figEJ 



#3 We want to thank M.R. Pennington for a detailed comparison of his results with ours and interesting discussions. 
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4 Conclusions 



We have undertaken a dispersive study of the 77 — > 7r°7T reaction. Our approach is based on that of 
refs. |18t \19\ [T7] but using a better behaved Omnes function for the 1 = S-wave irir channel. As a 
result, we have been able to reduce drastically the uncertainty regarding the 4>o(s) used in this Omnes 
function above the KK threshold. Our improvement is equivalent to take three subtractions instead of 
the two originally proposed in refs. |18l [T9l 117] . We have then used two low energy conditions and a third 
constraint in the form of a bound on the /o(980) region so as to fix the three subtraction constants. This 
has allowed us to present more accurate results, which might be used to discriminate between different 
7T7T 1 = S-wave parametrizations, once more precise data on 17(77 — * 7r°7r°) become available. Further 
improvements at the theoretical level rest on a more precise determination of 4>o( s ) above sk and a more 
systematic calculation of Lj(s). 

We have calculated the ratio of the residua \g all / gamr\ = (2.1 ±0.25) X 10~ 3 with s CT from refs. [45, 46] . 
The a width to 77 was also studied and we stressed its dependence on the s a and IsWt,-! 2 employed, not 
used in our dispersive study of 77 — ► 7r°7r°. One value obtained is T(a — ► 77) = (1.8 ± 0.4) KeV with 
s a and IsWttI — 3 GeV from ref . |45j . The others values take s CT as given in ref. [46] with T(a — > 77) = 
(2.1 ± 0.3) KeV for \g a7T7r \ = 3.4 GeV, and T(a -» 77) = (3.0 ± 0.3) KeV for \g a7T7r \ = 3.9 GeV. The last 
two numbers for T cd (a — ► 77) tell us that the uncertainties in its calculation are still rather large and a 
further improvement requires to know precisely \g^ n \ from ref.[37j. 
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